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Abstract 



We first study how sterile neutrinos can fit the 5cr z/„ 



LSND anomaly: 2+2 solutions 
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are strongly disfavoured by solar and atmospheric data, while 3+1 solutions can still give 
a poor fit (for a specific range of oscillation parameters, to be tested by MiniBooNE). If 
MiniBooNE will see no — > z/g transitions, we will have a hint for CPT violation. Already 
now, unlike sterile neutrinos, CPT-violating neutrino masses can accomodate all safe and 
unsafe data. We study how much CPT must be conserved according to atmospheric and 
K2K data and list which CPT-violating signals could be discovered by forthcoming solar 
and long-baseline experiments. 



Oscillations between the three Standard Model (SM) neu- 
trinos are described by two independent squared neutrino 
mass differences, allowing to explain only two of the three 
neutrino anomalies (atmospheric solar |2j and LSND P|) 
as oscillations. A joint fit is not possible even if one trusts 
only the safest data from atmospheric, solar and reac- 
tor ^ neutrino experiments: the the up/down atmospheric 
asymmetries and a ~ 50% disappearance of solar v^- Most 
global fits of neutrino data drop the LSND anomaly be- 
cause the other ones are considered as more solid. In quan- 
titative terms, we have a Scr solar anomaly (although it can 
be reduced to 5cr by dropping solar model predictions), a 
14cr atmospheric anomaly and a 5(7 LSND anomalyt. The 
'number of standard deviations' is here naively computed 
as (Ax^)^/2 = (xIm - Xbcst)^^^ where xLst is the 
value corresponding to the best-fit oscillation, and xIm 
corresponds to massless SM neutrinos. 

*In the addendum at pages I9lll2l wc update our results includ- 
ing the first data from KamLAND and WMAP, which disfavour the 
CPT-violating and '3+1' solutions. 

tOn leave from dipartimcnto di Fisica dcirUnivcrsita di Pisa and 
INFN. 

t The Dp^ — > Pe LSND anomaly is presented as an evidence for 
a. II ^ e oscillation probability of (0.264 ± 0.081)% 0, that differs 
from zero only by slightly more than 3cr. However, from a table of the 
likelihood C, obtained from the LSND collaboration and computed 
on an event-by-event basis, we read 
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rather than ~ 10. 



A reanalysis of LSND data that chooses stronger cuts obtains Ax^ = 
47 (eq. (2.7) of |51). These large Ax^ mean that the LSND anomaly 
cannot be due to a statistical fluctuation. It is not clear which data 
really contain the LSND evidence. Apparently, some mark of oscil- 
lations that cannot be summarized by the number of Pe events is 
hidden in the full LSND data, maybe in the energy distribution. 



In section n we discuss how and how well oscillations 
with extra sterile neutrinos can fit the LSND anomaly |H] . 
In particular we study which one of the two different kind 
of four-neutrino spectra (3+1 or 2+2) is favoured by the 
present data, and by an eventual future confirmation of the 
LSND data. Taking into account the recent SNO result (2] 
an extra sterile neutrino can improve the situation only in 
the 3+1 scheme, and even this case does not allow to fully 
reconcile all data. 

This situation suggests to look for alternative interpre- 
tations of the LSND anomaly. One possibility is that either 
the atmospheric or solar or LSND anomaly is not due to 
oscillations. Various mechanisms (even unplausible ones) 
can fit the data as well as oscillations |Z| IHj ■ 

Using only oscillations, all data can be consistently fit- 
ted by the CPT-violating neutrino spectrum illustrated in 
fig. n This solution was proposed in [H] when the initial 
2.6ct LSND hint for Vf. ^Hl decreased down to 0.6fT, 

leaving an anomaly only in i>f^ El ■ Unlike sterile neu- 

trinos, this solution also satisfies (unsafe?) bounds from 
nucleosynthesis and SN1987A 0111]. Despite the lack of 
theoretical grounds, this speculation is interesting because 
can be tested soon. If CPT violation were the right an- 
swer, MiniBooNE jE] (the experiment designed to test 
LSND, looking for i'^ — s- fe) will not see the LSND oscil- 
lations; a — > i/g experiment is needed to directly test 
this possibility. If CPT is badly violated as in fig. ^ one 
generically expects detectable CPT-violating signals in at- 
mospheric and solar oscillations. In any case it remains 
interesting to constrain CPT-violation in neutrino masses. 
In section |21 wc compute the present bounds and list the 
possible CPT-violating signals and surprises that could ap- 
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Figure 1: The CPT-violating spectrum proposed in )^. 



pear in forthcoming solar and long-basehne experiments. 



1 Sterile neutrinos 

The sterile neutrino can be used to generate either the 
LSND or the solar or the atmospheric anomaly. 

3+1 neutrinos 

Within this scheme the sterile neutrino is employed to gen- 
erate the LSND anomaly. In fact, in the jargon 3+1 in- 
dicates that the additional sterile neutrino is separated 
by the large LSND mass gap from the 3 active neutrinos, 
separated among them only by the small solar and atmo- 
spheric mass differences. A theoretical remark is in order. 
If the 4x4 neutrino mass matrix mai (i — {i, s} and 
£ = {e, /i, t}) has the naive form 



mis 



Toss = ?7lLSND, 



^ TOs 



the sterile neutrino induces a contribution to the solar mass 
splitting of order-'- 



^"^LSND 20LSND 



10-(3-i) eV^ 



that is too large in most of the region allo-wed by solar 
and LSND data. One needs either a cancellation or a 
mass matrix of the special 'approximatively rank one' form 

"T-ii' — ^is^i's'TlLSND- 

Even ignoring this potential theoretical problem, 3+1 
oscillations present a phenomenological problem, because 
predict that — > Ve oscillations at the LSND frequency 



proceed trough 



and 



are strongly 



constrained by disappearance experiments. More precisely, 
keeping only oscillations at the dominant LSND frequency 

S' = sin2(ATO^SNDi/4£^^) 



•fMore precisely, assuming 6ig <^ 1, maximal atmospheric mixing 
and 6i3 = 0, and taking into account the larger atmospheric mass 
splitting, one has SAm^^^ = {9^^ + 6^^)-^ Am^gj^j^ where 9j_i, 



S sin^ 29^ 
S sin^ 29^ 



P{v,^v,) = 1 
P{l^^.^^^.) = 1 
P{ve J^m) = ^ sin^ 26'lsnd 

with 6'lsnd ~ Oesdfj.s, or more precisely ^31 



sin^ 29 



LSND = ^ sin 26'esSin 26'^s- 



(1) 



The 9es mixing angle is constrained by Bugey, Chooz 0], 
SuperKamiokande (SK) atmospheric data P and the 9^s 
mixing angle by SK, CDHS and CCFR Furthermore 
— » i/e oscillations are also directly constrained by Kar- 
men ^H]. Fig. |21 illustrates ho-w accurately we reproduce 
such bounds^. 

The crucial question is if these bounds are too strong 
for allo-wing the oscillations suggested by LSND. At first 
sight the ans-wer is that they are ^3], but this negative 
conclusion -was questioned in |17 | and the first accurate 
statistical analysis of this issue -was performed in JS| with 
Bayesian techniques. Our result, shown in fig. |2| basically 
agrees with |18j . Working in gaussian approximation^ we 
find that all 96% CL LSND confidence region is excluded 
at, at least, 96% CL level. Therefore 3+1 solutions have 
some goodness-of-fit problem. One needs to invoke a sta- 
tistical fluctuation with around % probability to explain 
why only LSND sees the sterile oscillations. 

Even if this conclusion is self-evident, we justify the 
adopted statistical strategy. As discussed in due to 
the large number of d.o.f. (about 200) a naive Pearson 
global test is unable to notice this problem and would 
erroneously suggest that 3+1 oscillations give a good fit. 
While it is difficult to develop a general and efficient good- 
ness-of-fit test, in this particular case the fit is bad for one 
specific reason: different sets of data are mutually exclusive 
(up to a 96% CL) within our theoretical assumptions. In 
such a situation the goodness-of-fit problem is efficiently 

§Wc used the SK atmospheric results Q after 79 kton-year (55 
data), K2K 1161 (at the moment K2K finds 44 events, versus an ex- 
pected no-oscillation signal of 64 it 6 events), the latest solar results 
from Homestake, Gallex, SAGE, GNO, SK, SNO (49 data), the fi- 
nal BuGEY (60 data), Chooz (14 data), CDHS (15 data), CCFR (15 
data), Karmen and LSND results. We use the likelihoods computed 
by the Karmen and LSND collaborations on an event-by-event basis. 
We have not included data from MACRO Q (that confirms the atmo- 
spheric anomaly) and from earlier atmospheric experiments because 
are less statistically significant than SK. The data are combined by 
multiplying all likelihoods C (i.e. by summing all = — 21n£). At 
Am^ ^ 10 eV^ Chooz and BuGEY bounds could bo considered as not 
fully trustable because limited by the theoretical error on the total 
Pe fluxes generated by reactors. 

^So that Ax^ = 7 corresponds to 97% CL level for the two pa- 
rameters ^^LSND S'lid ^"^LSND' '^'^^ Gaussian approximation is not 
fully satisfied (e.g. our best fit regions are not ellipses). A Bayesian 
analysis can shift 97% to ~ 95% or ~ 98%, with 'reasonable' choices 
of the prior probability distribution, (the arbitrarily remains until 
there are 'large' allowed regions). As discussed in 1191 . a similar shift 
is typically obtained in a frequcntist analysis, that cannot however 
be performed in a reasonable computing time. Therefore we stick to 
the Gaussian approximation. 
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Figure 2: 90% CL regions from Karmen, CDHS, CCFR, 
BuGEY, Chooz and LSND (shaded). The mixing angle 6 
on the horizontal axis is different for the different experi- 
ments. 



Figure 3: The LSND region at 90% and 99% CL, compared 
with the 90% (dashed line) and 99% CL (continuous line) 
combined exclusion bounds from data in fig. [3 and SK. 



recognized by fitting separately the two incompatible data. 
This is what is done in fig. |31 

Ignoring the poor quality of the fit, the best combined 
fit region for the LSND parameters is shown in fig. It 
agrees reasonably well with the corresponding fig. in |20j . 
taking into account that we show values of 

X^(^LSND, AmLSNo) = minx^b, 6'lsnd, Ajtzlsnd) 
p 

(where p are all other parameters in which we are not 
interested), so that we convert values of — Xbost i'^to 
confidence levels using the gaussian values appropriate for 
2 d.o.f. (the 2 LSND parameters), while a statistically less 
efficient procedure with more d.o.f. is employed in |20|. 

2+2 neutrinos 

In the jargon 2+2 indicates 2 couples of neutrinos (one 
generates the solar anomaly, and the other one the atmo- 
spheric anomaly), separated by the large LSND mass gap. 
Within this scheme, the sterile neutrino is employed to gen- 
erate the solar or atmospheric anomaly, or one combina- 
tion of the two. The fraction of sterile neutrino involved in 
solar oscillations, r?!"", plus the fraction of sterile neutrino 
involved in atmospheric oscillations, rjf"™, is predicted to 
sum to unity ^T) 

Experiments now tell that both the solar and atmospheric 
anomalies are mostly generated by active neutrinos, and 
only a small sterile contribution is allowed. Consequently 
2+2 oscillations give a global fit worse than 3+1 oscilla- 
tions 122 ED]. Let us summarize the present experimental 
status of this issue. 



• Solar data give a 5.4(t evidence for pure active solar 
oscillations versus pure sterile oscillations: combin- 
ing all solar data in a global fit we obtain " 

Xsun(best sterile) — Xsunlbcst active) — 30 

and 77f " = ± 0.18. In particular, SNO/SK find a 
5.1(7 direct indication for z^^ appearance. 

• Atmospheric data data give a la indication for 
pure active atmospheric oscillations versus pure ster- 
ile oscillations. In fact, a global fit of atmospheric 
data gives [111^]** 

Xatm(best sterile) - Xatm(best active) « 50 

and rif"™ = 0+0.16. This strong evidence is obtained 
combining independent sets of data. SK claims 
that pure sterile is disfavoured by the up/down ratio 
in a NC-enriched sample (3.4 standard deviations) 
and by matter effects in partially contained events 

II Some words of caution. Arbitrary choices become more relevant 
when fitting disfavoured data (for example: the error is evaluated at 
the experimental point or at the theoretical point?). Furthermore, 
our bound on the sterile fraction allowed by solar data is obtained 
assuming the BPOO 1231 prediction for the Boron Ve solar flux. It is 
proportional to the '^Bep — > *B7 cross section: some authors think 
that systematic uncertainties in its measurement could be underes- 
timated. 

**A large amount of these atmospheric data is not included in 
theoretical reanalyses (because not yet accessible outside the SK col- 
laboration in a form that allows to recompute them) that therefore 
obtain a much smaller Ax^ 15 HOH^ in place of 50 This 
underestimation of the SK bound means that at the moment only 
SK can perform a sensible analysis of mixed sterile and active at- 
mospheric oscillations and explains why the authors of I2UI do not 
recognize that 2+2 oscillations are extremely disfavoured. One mix- 
ing angle is set to zero in the SK analysis; relaxing this unjustified 
simplification should not significantly weaken the bounds. 



3 



10^ 



10 



10" 




LSND as Vm 



10" 



10^ 



10 



10" 



10" 




10" 



10"-^ 



10" 



10" 



10" 



10" 



sin^lOLs 



10"" 
sin^2eLSND 



10" 



Figure 4: Best-fit regions at 90% and 99% CL (2 d.o.f.) for the LSND parameters assuming oscillations. 

assumes that the LSND anomaly is generated trough a sterile neutrino ('3+1' scheme). Fig.^) assumes that 
the LSND anomaly is generated by active neutrinos, while something else (e.g. neutrino decay, sterile neutrinos,. .. ) 
generates either the atmospheric or the solar anomaly, without affecting LSND. The dotted lines show the regions 
suggested by only the LSND data. The dots show the best fit points. 



2.9a) and upward through-going muons (« 2.9(t). 
In total 7(7 P EI]. Matter effects in MACRO [T] give 
another 3.1a signal. Furthermore SK finds a direct 
2a hint for r appearance. 

In summary, the two extreme cases (all the sterile in atmo- 
spheric oscillations and all the sterile in solar oscillations) 
have been excluded, as summarized in table At the mo- 
ment published results only allow an approximated analy- 
sis of intermediate cases. We find that 77^°* = 0±0.25, with 
j^tot _ ^ disfavoured at 4a. Intermediate cases are less dis- 
favoured than the two extreme cases by only the amount 
expected, on a statistical basis, due to the presence of one 
more parameter: the 'best' fit is now obtained around the 
weighted average of the two incompatible solar and atmo- 
spheric determinations, ry^*™ = 1 — « 0.5. We do not 
present more precise results because fitting incompatible 
data makes little sense. Despite the approximation, the 
final conclusion is clear: 2+2 oscillations are too strongly 
disfavoured to be considered as a viable possibly. 

In fig. 2Jd we show the best-fit region for the LSND 
parameters, assuming that the LSND anomaly is generated 
by oscillations of active neutrinos. This result applies to 
a general class of models where something different than 
oscillations between active neutrinos is the source of the 
solar or atmospheric anomalies. In particular it applies to 
2-1-2 oscillations: despite they are strongly disfavoured the 
LSND best-fit regions are unaffected by the problems with 
solar and atmospheric data, and can therefore be reliably 
computed. 

This region extends to values of the LSND parameters 
not accessible within 3-fl oscillations, see fig.E] Therefore 
the value of — > i^e) that will be measured at Mini- 



BooNE could discriminate between the two cases: roughly, 
3-1-1 oscillations prefer a value of P{Vf_i Ve) somewhat 
smaller than the one suggested by LSND. Furthermore 3-1-1 
spectra must be accompanied by a significant disappear- 
ance of at the LSND frequency. For example, our 3-1-1 

best-fit (marked with a dot in fig. Ela,) has sin^ 29 ~ 0.2, 
around the sensitivity of MiniBooNE. 

Both 2-1-2 and 3-1-1 oscillation patterns can be real- 
ized with different neutrino spectra. Since at the moment 
(and in the near future) no experiment can resolve the dif- 
ference we do not consider all possibilities. For example, 
even knowing the oscillation parameters and the type of 
spectrum, we could not safely predict neutrinolcss double 
/3 decay signals. 



Many sterile neutrinos 

As shown in the last paper in |17| . many sterile neutrinos 
cannot give a much better 3-1-1 fit than a single sterile 
neutrino. Of particular interest are minimal models where 
right-handed neutrinos live in a single extra dimension of 
radius R [^Hl, that could be identified with the LSND scale. 
In such 3 -|- cxD models the problematic prediction of 
3+1 oscillations becomes slightly more problematic |Hj. In 
fact, for small mixing angles and in the limit of averaged 
sterile oscillations, we now have ^lsnd ~ a/T/IO 0es(^fj.s in 
place of 0LSND ~ Oes(^fj,s- More importantly, the effective 
active/sterile mixing angles are now predicted to be 



= Y^Val'Aml^R^ 
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(for a hierarchical spectrum of active neutrinos, the other 
cases are more problematic). The Chooz bound on Ve3 
(that will soon be tested and eventually strengthened by 
long-baseline experiments) now gives another constraint 
on 6es , making this minimal model more problematic than 
3+1 oscillations. One can consider a large variety of less 
predictive non-minimal extra dimensional models. 

In the case of sterile solar or atmospheric oscillations, 
many sterile neutrinos can be less disfavoured that a sin- 
gle sterile neutrino. As discussed above, pure atmospheric 
sterile oscillations are disfavoured mostly by matter effects 
(in the earth), that suppress i'^ — > Vg at large energy: 
SK data are better fitted by z/^ — > i/^ oscillations, unsup- 
pressed by matter effects. Even in the solar case, mat- 
ter effects (in the sun) contribute to determine how much 
SMA sterile oscillations are disfavoured [201 • In presence 
of a tower of many sterile neutrinos, matter effects do not 
suppress sterile oscillations at large energy or density, un- 
til there is a sufficiently heavy sterile resonance to cross. 
However, sterile oscillations must be strongly matter sup- 
pressed within a supernova. As discussed in |S] super- 
novae strongly constrain sterile towers that continue up 
to masses of 10''^"^^ eV. This is e.g. the case of an extra- 
dimensional Kaluza-Klein tower that continues up to the 
TeV scale In conclusion, (2 + many) oscillations can 
be less disfavoured than 2 + 2 oscillations. However, even 
forgetting the lack of theoretical motivation, it docs not 
seem possible to achieve a really satisfactory fit. 

2 CPT violation 

Theory 

The only safe result is that CPT is conserved in Lorentz- 
invariant local quantum field theories (QFT). Therefore 
CPT-violating effects can be obtained by abandoning lo- 
cality or Lorentz invariance: 

1. In local QFT, CPT violation can be induced if the 
Lorentz symmetry is broken, e.g. spontaneously by 
vacuum expectation values of fields with spin 1 or 
higher, or cosmologically by interactions with some 
'aether', or by a non-trivial extra-dimensional back- 
ground, or... 

This first possibility |27II28| seems not promising for LSND: 
like anomalous matter effects and unlike oscillations, new 
effects arc not enhanced at low neutrino energy. There- 
fore old experiments ^2 done at energies 2 -f- 3 orders of 
magnitude higher than LSND, disfavour the best fit Kar- 
men/LSND region. Furthermore in this context it seems 
difficult to obtain Pee < 1/2 (as suggested by the latest 
SNO data E]) in solar oscillations |30| . 

Therefore we focus on the second possibility, that could 
explain the LSND anomaly P|: 

2. Strings, branes, quantum foams, wormholes, non com- 
mutative geometry (and other non local things like 



that) suggest CPT-violating effects, maybe suppres- 
sed by only one power of the quantum gravity scale 
M (this case gives rise to interesting signals even for 
M - 10^9 GeV ED)- 

If an effect at that level were an unavoidable phenomenon, 
quantum gravity at the TcV scale would be excluded by 
bounds on the KqKq mass difference: 

niKo -^Ka < 0-4 10~^eV. 

The mass difference between neutrinos and anti-neutrinos 
that could explain LSND is larger by many orders of mag- 
nitude: we assume that CPT-violating effects are domi- 
nantly felt by neutrinos. 

The generic Hamiltonian that describes non relativis- 
tic systems (e.g. Kaons) violates CPT, if the constraints 
from the underlying local relativistic QFT are not im- 
posed. In the case of relativistic systems (e.g. neutrinos) 
one can mimic the standard Hamiltonian demanded by lo- 
cal relativistic QFT (particles together with anti-particles) 
but without imposing all the constraints demanded by 
QFT (particles degenerate with anti-particles), so that the 
generic Hamiltonian that describes free propagation of Dirac 
neutrinos has different mass terms for v and P. The social 
duty of studying how CPT-violating neutrino masses can 
arise in popular fundamental models has been exploited 
in |32j ■ obtaining the imprimatur from string branc- world 
orbifolds. Non commutative geometry was invoked in Pj. 

We do not consider other possible CPT violations in 
neutrino interactions, because experiments with (mainly) 
i^^ , Pfj, beams and precision electroweak data [HSl 1^ find 
that neutrino NC couplings cannot differ from the SM pre- 
diction by more than few %. A global fit of electroweak 
precision data |34| shows that the CC couplings of e and /i 
neutrinos agree with the SM with few per-mille accuracy. 

Fit of SK and K2K data 

In absence of oscillations, the number of I'^-induced events 
at SK would be roughly double than the number of i/^- 
induced events (the ratio is higher at sub-GeV energies. 
This is mainly due to the different and 9^^ cross-sections 
on matter, that we compute by summing the elastic and 
deep-inelastic cross sections [23 )■ We assume that SK has 
an equal efficiency for v and i/-induced events. 

We use a (hopefully) self-explanatory notation for the 
v and D parameters. An over-bar marks anti-neutrino pa- 
rameters. For example, ^atm and Afh'^^^ parameterize the 
atmospheric — > Dt oscillations. 

Restricted analysis To begin, we assume that 6*8^1, 
dcKOOZ, ^CHOOZ, ^LSND havc negligible effect on atmo- 
spheric oscillations^ that are therefore described by Am^^j^, 
Am^j^, 6'atm and 0atm- 

A simple approximation captures the main properties 
of the fit. The up/down asymmetry in the number of 
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multi-GcV muon events is pP 

N, - N't 
A = ^ = 0.327 ± 0.045 

Assuming maximal mixings, in the CPT-conserving case 
one has 



-2 ~3 10-^eV^: Awl/3 (2a) 
The asymmetry is smaller in CPT-violating cases, e.g. 



Am: 



atm ^ ^'^atm 



3 10" 
3 10" 



'eV^ 
'eV^ 



3 lO^^eV^ 
3 10"^ eV^ 



A«l/4 (2b) 
Awl/5 (2c) 
(2d) 

A « 1/11 (2e) 



and even smaller if mixings are non maximal. These con- 
siderations allow to understand the main features of our 
numerical result. In fig. [3 we show the minimized with 
respect to the mixing angles 6'atm and ^atm- While Am^^j^j 
is almost as strongly constrained as in a CPT-conserving 
fit, Am^jj^j can be about one order of magnitude larger or 
smaller that Amatm-^^ The global for SK data is here 
obtained by summing the corresponding to the individ- 
ual zenith-angle distributions of sub-GcV and multi-GeV 
(10 e-like bins and 10 fi-likc bins each), stopping fi (5 bins) 
and upward-through-going /x (10 bins) events. The overall 
normalization in each kind of events has been considered 
as a free parameter. 

Alternatively, one can try to take into account the theo- 
retical predictions for the overall fluxes as in j38j employing 
a 55 X 55 correlation matrix. This second approach gives 
a slightly different bound on CPT-violation: larger values 
of 

^™atm would not be significantly disfavoured up to the 
right border of fig. El 

Since the best fit is obtained for almost CPT-conserving 
oscillations, the fit for the mixing angles is quite simple, 
and we do not need to show a dedicated figure. In the 
CPT-conserving case sin^ 20atm has to be close to one. We 
find that in the CPT-violating case the same bound applies 
replacing 



sin'' 261a 



sin^ 20atm 



sin^ 20atm 



so that both O^tm and (to a lesser extent) 6'atm have to be 
close to maximal. 



ttAn analogous fit of sub- and multi-GeV SK data has been per- 
formed in I.Sfil . finding Ax^ = Xcpt~Xc*^ = 16, while we do not find 
any strong evidence for CPT-violation. As clearly discussed in 1,361 
this large Ax^ could be an artifact due to having neglected the error 
on the ratio between and Ue fluxes. Our results also disagree with 
another CPT-violating fit presented in 1281 : the difference is signif- 
icant even in the CPT-conserving limit. A fit performed by the SK 
collaboration 1211 agrees with our fig. |21 

In the case of K2K data (sensitive to neutrinos) we fitted the total 
number of events ignoring the information about their energy, finding 
a result in agreement with 1371 . 



General analysis We now discuss the effects of the other 
mixing angles, 6sun, ^'chooz, ^chooz, ^lsnd, that we have 
so far neglected. Some of them are allowed to be large, but 
cannot significantly affect our CPT-violating atmospheric 
fit shown in fig. 

In anti-neutrinos, disappearance experiments require 
small values of the two mixing angles that induce oscil- 
lations at the LSND frequency. These constraints allow 
for a novel possibility, somewhat disfavoured only by at- 
mospheric data: the most splitted anti-neutrino eigenstate 
could be dominantly i/^ (rather than i?e as in fig.Ul). In this 
case, ^cHOOZ (the remaining mixing angle that now gives 
oscillations at the atmospheric frequency) could be large, 
without confiicting with the Chooz bound, if Arhl^^-^ is 
below the Chooz sensitivity. 

In neutrinos, solar experiments require ^sun ~ 1 as in 
the CPT-conserving case. Unlike in the CPT-conserving 
case Chooz does not force Am'^^^<0.7 10~^eV^, but a 
larger Am^^^ has recently been disfavoured by the SNO 
NC data The angle ^chooz (that induces i/^ Vf. 
oscillations at the atmospheric frequency; we improperly 
adopt the name used in CPT-conserving analyses) is not 
bounded by Chooz (i.e. by disappearance of Ve)-, but only 
by global fits of solar and atmospheric data, that weakly 
disfavour a large ^chooz 

Signals At the light of these results, we can now list the 
CPT-violating signals that could appear in forthcoming 
experiments (some signals were discussed in 1321 



• MiniBooNE will not see the LSND oscillations, if will 
only search them as v^^ 



Vf. rather than as t'p — > 



While this signal is mandatory if the CPT-violating inter- 
pretation of the LSND anomaly is correct, the following 
signals can but need not to appear, depending on the val- 
ues of the unknown parameters: 

• We would have a signal for CPT violation if Kam- 
LAND will find no solar oscillations in its reactor 
data, and Borexino will indirectly favour LMA by 
finding a ~ 1/2 suppression and no matter nor sea- 
sonal effects. 

• If ^CHOOZ were large, KamLAND would discover 
its effects and misinterpret them as LMA oscilla- 
tions. In particular this means that if KamLAND 
will confirm LMA, a CPT-violating interpretation of 
the LSND anomaly would not be immediately ex- 
cluded, but only disfavoured. We do not list other 
possible situations that could happen depending on 
future Borexino and KamLAND results. M 

• According to our fit in fig. |S1 long-baseline exper- 
iments that plan to employ a i/^ beam (like K2K, 



"'"'"When these results will be announced, we will update the hep-ph 
version of this paper, adding a precise discussion. 
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Table 1: Interpretations of all oscillation data, ordered according to the quality of their global fit. The last 
4. columns show the minimal Ax^ restricted to solar data, atmospheric data, LSND data, and to experiments compatible 
with no oscillations (mainly Chooz, Bugey and CDHS). 



Minos and CNGS) have almost the same capabili- 
ties of confirming atmospheric oscillations as in the 
CPT-conserving case. Using a. a. beam they can 
also test if Am^j^ is higher than Am^^,,^ (if Ajfi^jj^^ 
is as large as possible, a 5% contamination in the 
beam could also give detectable r-appearance ef- 
fects). 

• These long-baseline experiments can test if ^chooz 
is larger than what allowed in the CPT-conserving 
case by looking at ~* Ve- 
in longer terms, an atmospheric experiment that sepa- 
rately measures Arn^^^ and Am^^j^^ (and sees the first os- 
cillation dip) seems feasible 00], although KEK, CERN 
and FermiLab preferred to pursue 3 long-baseline experi- 
ments. 

With a hierarchical v spectrum (rather than with the 
inverted spectrum motivated in [Hj) planned /3-decay ex- 
periments like KATRIN 0Tj can test the upper part of the 
Am^ range suggested by LSND Planned neutrino- 

less double /3-decay experiments 021 have brighter perspec- 
tives of improvement than /3-decay experiments, but CPT- 
violating neutrino masses seem to require Dirac (rather 
than Majorana) neutrinos, if the Lorentz symmetry is un- 
broken (because there is no Lorentz-invariant distinction 
between massive Majorana v into a, D: a sufficiently 'fast' 
Lorentz transformation transforms v in P). 

In the far future, with a neutrino factory it should be 
possible to test CPT conservation in atmospheric oscilla- 
tions at the % level |43j . 

3 Conclusions 

A possible global explanation of the three neutrino anoma- 
lies (atmospheric, solar and LSND) is that an extra sterile 
neutrino generates one of them. Each anomaly, when fitted 
independently from the other ones, prefers active oscilla- 
tions refusing the sterile neutrino. The relatively better 
global fit is obtained with a 3+1 spectrum (sterile LSND 
oscillations) rather than with a 2-1-2 spectrum (sterile solar 
or atmospheric oscillations: this case is disfavoured at 4cr, 
after the recent SNO NC results |21). However the fit is 
not good: within the 3-f 1 scheme the LSND anomaly con- 
flicts with z^e or Vfj^ disappearance experiments. One needs 



to invoke a statistical fluctuation with around % probabil- 
ity to understand why Bugey, Chooz, CDHS or SK have 
not seen sterile effects. Our main results are summarized 
in table 1. 

The best-flt LSND regions are shown in flg. 0] assum- 
ing that the LSND anomaly is generated trough a sterile 
neutrino (3-1-1 case, flg.EK) or by oscillations of active neu- 
trinos (flg. ^p) , assuming that a sterile neutrino or some- 
thing else (e.g. neutrino decay) generates the solar or atmo- 
spheric anomaly. The best flt LSND regions are somewhat 
different: MiniBooNE could discriminate the two cases. 

Many sterile neutrinos (motivated e.g. in extra dimen- 
sional models) can somewhat improve the fit, but it does 
not seem possible to obtain a good sterile solution. 

In view of these unsatisfactory sterile fits, and of the 
latest LSND results 

Piiy^^i^e) = (1.0 ±1.6) 10-3 
Piv^^^ve) = (2.6 ±0.8) 10-3 

one might want to speculate on CPT-violation. A sat- 
isfactory global fit of all neutrino data (see table 1) can 
be obtained with the CPT-violating neutrino masses pro- 
posed in in]. Theory gives no useful restriction, and in 
particular does not tell if CPT should be violated also in 
atmospheric oscillations, although it looks plausible. Fig.|2] 
shows how present SK and K2K data restrict the atmo- 
spheric oscillation parameters Am^j^^ and Affi^^^. They 
can differ by about one order of magnitude. In section 2 
we studied which CPT-violating oscillations are compati- 
ble with present data, and listed the unusual signals that 
could be seen at forthcoming solar (KamLAND, Borexino) 
and long-baseline experiments (K2K, MINOS, CNGS) — 
and of course at MiniBooNE. 

Acknov^rledgments We thank K. Eitel and M. Steidl 
for giving us the table of the Karmen likelihood, and B. 
Louis and G. Mills for the LSND hkelihood. We thank V. 
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Figure 5: Fit of SK and K2K data for the neutrino 
and anti-neutrino atmospheric mass splitting at 68, 90, 
99% CL (2 d.o.f). 

4 Addendum about the first Kam- 
LAND and WMAP data 

CPT-violating solution The CPT-violating neutrino 
spectrum suggested in [SJ allowed to reconcile the so- 
lar, atmospheric and LSND neutrino anomalies but pre- 
dicted no effect in KamLAND. The evidence seen by Kam- 
LAND [l^ can be fitted by the alternative CPT-violating 
spectrum proposed at page of the original version of this 
paper, at the price of a non standard fit of atmospheric 
data. We now show that, as anticipated in j2 A| . this solu- 
tion is disfavoured by atmospheric data (we disagree with 
the claim in |3Aj that it gives an atmospheric fit "clearly 
favored over the CPT conserving one"). 

We denote anti-neutrino parameters with an over-bar. 
Fitting SK and K2K data in the usual two neutrino ap- 
proximation, in fig.Elwe plotted the best-fit regions for the 
atmospheric mass splittings (Am^^j^, Arri^^^), marginal- 
izing the global with respect to the atmospheric mix- 
ing angles ^atm and 0atm- We can now do the same test 
on solar u data and reactor D data (from the Homestake, 
SAGE, Gallex, GNO, SK, SNO, KamLAND, CHOOZ ex- 
periments). The result is shown in fig. El The plot is re- 
stricted to the LMA solution: other solutions with smaller 



Figure 6: Fit of solar and reactor data for the neutrino 
and anti-neutrino solar mass splitting at 68, 90, 99% 
CL (2 d.o.f.). 

ATOg^jj (LOW and QVO solutions) are disfavoured but 
not excluded by solar v data. The best fit is close to 
CPT-conscrvation in both the atmospheric and in the so- 
lar cases. We do not show the corresponding fits for the 
mixing angles 6'atm, ^atm, 6'sun and O^nn- AH these mixing 
angles must be large. 

After including KamLAND data, CPT-violation can no 
longer perfectly fit all data. We now have some evidence 
of 'solar' and 'atmospheric' oscillations not only in v but 
also in v, leaving no room for the larger Am^ that should 
give rise to the LSND anomaly. 

In both 'solar' and 'atmospheric' cases v data do not 
yet provide a conclusive evidence, so that we may attempt 
to fit all neutrino anomalies by sacrificing either i) solar 9 
data or ii) atmospheric u data. 

Concerning case i), we just mention that the CPT- 
violating spectrum proposed in jU] I32| (that predicted no 
anomaly in KamLAND) still gives a reasonably good global 
fit. We do not consider this possibility, that KamLAND 
should exclude with more statistics. 

Rather, we explore case ii) and sacrifice atmospheric 
data, taking the smaller anti-neutrino Am^ in the Kam- 
LAND (KL) range, rather than in the atmospheric range, 
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and the larger Ato^ in the LSND range. There are three 
anti- neutrino mixing angles. Two of them (that from now 
on we name ^lsnd and ^atm) give oscillations at the larger 
LSND frequency, Atolsnd- ^he third mixing angle (that 
from now on we name ^kl) gives oscillations at the smaller 
Am^L. KamLAND and LSND data want a large Okl ~ 1 
and a small ^lsnd ~ (0.1 ^0.01). The last mixing angle, 
^atm, induces <-> 9^ oscillations at the LSND frequency. 
As usual, the relatively better atmospheric fit is obtained 
for maximal ^atm- 

Beyond performing a global fit, it is useful to present 
a semi-quantitative understanding of SK data. For our 
purposes the main observables are the number of up-ward 
going and down- ward going /i-like events in the multi-GeV 
sample. We recall that SK cannot distinguish ly from 9. 
Neutrinos have roughly the same flux and a two times 
larger cross section than anti-neutrinos: Analytical esti- 
mates are performed by just using this factor two. The 
total flux has a large overall uncertainty. 

We first give an argument similar to the one used in |3A| 
to state that atmospheric data favour CPT-violation, so 
that a comparison shows the reason of the disagreement. 

1. The maximal up/down asymmetry that CPT-conser- 
ving i/^ v-y oscillations can produce in the 'multi- 
GeV /i-like -I- PC sample of SK atmospheric data is 
(eq. ISEl)) 



^idcal 



jV^(cost9 = 1) - N^,{co'&^ = -1) 
7V^(cos?9 = 1) +iV^(cosi9 = -1) 



1 



where d is the zenith angle and cos = 1 corresponds 
to vertical down-going events. The corresponding 
upward/downward asymmetry is 

_ iV^(cosi9 > 0.2) - iV^(cosi9 < -0.2) 



A 



real 



0.28, 



7V^(cosi? > 0.2) -f iV^(cosi? < -0.2) 

defined ignoring 'horizontal' events with |cos6'| < 
0.2. The most recent SK data |4Aj give the value 



A. 



real 



0.288 ±0.030. 



The proposed CPT-violating scenario can give at 
most an up/down asymmetry Ajdcai = 1/4, which 
corresponds to Aroai = 0.21, l.ha below the experi- 
mental value. 

This argument takes into account only a part of SK data; 
furthermore the precise value of Ajoai /^^idoai depends on 



Ato^jjjj. We present another argument which avoids these 
drawbacks: 

2. Since Am^gND ^ ^"^atm ^^^^ proposed CPT-violating 
oscillations give a reduction in the muon rate with an 
energy and zenith-angle dependence which (up to an 
overall factor that does not play an important role in 
the SK analysis) can be mimicked by normal CPT- 
conserving oscillations with an appropriate effective 
value of the mixing angle 



sin^ 26*3 



^atm 
CPT — conserving 



4sin^ 26>at, 
6 - sin^ 2^a 



< 



Ami? 



The maximal value of the effective CPT-conserving 
sin^ 20atm allowed by the proposed CPT-violating os- 
cillations is 4/5, which is 5ct below the experimental 
value 

sin^ 26'at,n- = 1.00 ±0.04 Wt 



Our global fit of SK data, shown in fig. [SJ gives 

confirming the second argument. We remark that we fit 
all SK data and not only the up/down asymmetries. We 
computed the Ax^ using the latest SK data (1489 days of 
data taking''' |4A|'). Our numerical code contains precise 
neutrino and anti- neutrino fluxes and cross-sections. 

L-violating muon decay We extend our analysis con- 
sidering one more tentative interpretation of the LSND 
anomaly, in terms of a speculative AL = 2 muon decay 
channel with branching ratio roughly equal to the oscilla- 
tion probability suggested by LSND: 



(2.6 ±0.8) 10" 



Since the fact that LSND has a longer path-length than 
Karmen plays no role according to this interpretation, one 
naively expects that it is disfavoured by Karmen as much 
as oscillations with large Am? . 



This expectation was questioned in |5A| , that presented 
one explicit model that produces the /I — > €.9^9 decay with 

tWc thank M. Shiozawa for providing us the data presented at the 
Neutrino 2002 conference. Using the sHghtly older SK data set em- 
ployed in the original version of this paper would make no significant 
difference. 
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Figure 7: Status of 3+1 oscillations including WMAP data. Fig. EJa; LSND favours the shaded region. Values 
of AmLgND O'bove the horizontal dashed line are disfavoured by WMAP. The other dashed line shows the upper bound 
on 9 from all other neutrino experiments. The continuous line shows the combination of the two previous constraints. 
All bounds are at 99% CL for 2 dof. In fig. 0* we show the best fit 3+1 solution, including all data. 



Michel parameter p = (while p = 3/4 in ordinary muon 
decay) and consequently a Ue spectrum softer than the De 
spectrum produced by oscillations with large Am^. As a 
consequence the Karmen bound on appearance gets 
relaxed \5A\ by a factor A = 1.9 |6Aj -^ with respect to the 
bound obtained from the analysis in terms of oscillations 
with large Am^. In fact, Karmen detects using the 
DeP — > e?^ reaction, which cross-section is roughly propor- 
tional to E^^ . 

However also the LSND experiment detects De using 
the DeP — * en reaction. Therefore an interpretation of 
the LSND anomaly needs a BR(/i ev^v) larger than 
P{Vf^ Pg) by the same factor A. The two correction 
factors compensate each other when comparing LSND with 
Karmen, indicating that the naive expectation is right. 
Table |21 quantifies how much the /2 — » ev^v solution is 
disfavoured. A fully precise result would need a dedicated 
analysis, that only the LSND collaboration can perform. 

2+2 sterile neutrinos Rcf s |7AI 18 A| questioned the 
conclusion that an interpretation of the LSND anomaly 
in terms of an extra sterile neutrino with '2+2' spectrum 
have been excluded by solar and atmospheric experiments. 

■tWe thank Klaus Eitcl for communications about the Karmen 
analysis. 



The author of |7A| estimates that 2+2 oscillations pro- 
vide a global fit of all neutrino data with x^/dof « 291/276, 
which is acceptable. This is true, but the goodness-of- 
fit (gof) test based on the value of the total is inef- 
ficient when dof S> 1: it may assign an acceptable gof 
probability to a solution which is already excluded. This 
issue was discussed in the context of analyses of solar 
data in |2Aj and can be exemplified by recalling that, ac- 
cording to CPT-conserving global fits of solar and Kam- 
LAND data (e.g. |2Aj '). the LOW solution have been ex- 
cluded but its naive gof is still acceptable (presently it has 
xVdof = 89/91). 

The authors of |8A| pointed out that the full 4x4 mix- 
ing matrix contains small mixing angles which could relax 
the solar and atmospheric constraints on 2+2 oscillations. 
We do not expect that these effects can significantly im- 
prove the 2+2 status. 

3+1 sterile neutrinos After the first WMAP data cos- 
mology gives the dominant bound on neutrino masses |9Aj 

^771^ < 0.69 eV (95% CL, 1 dof). 

We assume that the extra sterile neutrino has a thermal 
abundancy (a possibility still compatible with primordial 
nucleosynthesis, unless uncertainties are aggressively esti- 
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model and number of free parameters 
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30 
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50 


SK atmospheric 



Table 2: Interpretations of all oscillation data, ordered according to the quality of their global fit. A 

= "f^ roughly signals an incompatibility at n standard deviations. 



mated jlOAp . so that its mass must be included in the 
sum. Both the WMAP bound and the assumed sterile 
abundancy are valid in standard cosmological model but 
might not hold in general. Present data are consistent with 
the standard cosmological model, and disfavour reasonable 
concrete mechanisms of evading the bound on J2 '^v (for 
recent discussions see e.g. [lOAIIllA] '). Therefore, keeping 
in mind the possible caveat that the bound on ^ mn could 
be weaker in some ad hoc cosmological model, we add this 
bound to our data-set. Extracting the precise probability 
distribution for ^ from |9A| , we study its consequences 
for LSND. 

As already remarked in jllAI ll2Aj (which performed 
semi-quantitative analyses and do not fully agree among 
them), the WMAP bound disfavours the 3-1-1 solution, 
since it employs a sterile neutrino with mass (AmLg^j))^/^ - 
^ m^. The precise combined bound is shown in fig.s0 ex- 
plained in their captions. These figures update our older 
fig.sOland^, and section^describes in greater detail how 
they are obtained. Table El quantifies how much the 3-1-1 
solution is now disfavoured. 

We do not study the impact on WMAP data on ster- 
ile solutions with 2-1-2 spectrum, since they are anyhow 
incompatible with neutrino data. 
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Conclusion We collect in table |21 the present status of 
various global interpretations of the solar, atmospheric and 
LSND neutrino anomalies. None of them allows to recon- 
cile all neutrino data in a clean way. It will be interesting 
to see if MiniBoone will confirm LSND. 
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